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We predict pp elastic differential cross section at LHC at the cm. energy y/s = 14 TeV and 
momentum transfer range \t\ =0-10 GeV 2 , which is planned to be measured by the TOTEM 
group. The field theory model underlying our phenomenological investigation describes the 
nucleon as a composite object with an outer cloud of quark-antiquark condensate, an inner 
core of topological baryonic charge, and a still smaller quark-bag of valence quarks. The model 
C****" , satisfactorily describes the asymptotic behavior of a to t(s) and p(s) as well as the measured 

pp elastic dcr/dt at -Js = 546 GeV, 630 GeV, and 1.8 TeV. The large |t| elastic amplitude 
of the model incorporates the BFKL Pomeron in next to leading order approximation, the 
perturbative dimensional counting behavior, and the confinement of valence quarks in a small 
■ region within the nucleon. 
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pp elastic differential cross section at LHC in near forward direction at cm. energy y/s = 14 
TeV and momentum transfer \t\ = - 10 GeV 2 is planned to be measured by the TOTEM 
(TOTal and Elastic Measurement) group [1]. Various models have been proposed to describe pp 
^ elastic scattering in the diffraction region \t\ ~ - 0.5 GeV 2 , such as: i) single Pomeron exchange 

with a trajectory ap(t)=1.08+0.25t [2], ii) multiple Pomeron exchanges with single- and double- 
diffractive dissociation[3], iii) the incident proton viewed as made-up of two color dipoles in the 
target proton rest frame [4]. pp elastic dcr/dt at LHC all the way from \t\ — to 10 GeV 2 has 
*^ been predicted on the basis of three different models: a) impact- picture model [5] based on the 

Cheng- Wu calculation of QED tower diagrams [6] , b) eikonalized Pomeron- Reggeon model using 
conventional Regge approach, but with multiple Pomeron- Reggeon exchanges included [7,8], c) 
effective field theory model that describes the nucleon as a chiral-bag with a quark-antiquark cloud 
[9,10]. A QCD motivated eikonalized model has also been proposed to predict pp dcr/dt at = 
14 TeV for \t\ = - 2.0 GeV 2 [11]. This wide array of models attempting to describe pp elastic 
scattering at LHC reflects the view that quantitative understanding of this process will provide 
fundamental insight into the nonperturbative and the perturbative QCD dynamics. 

The impact-picture model and the eikonalized Pomeron- Reggeon model predict besides the first 
dip-bump structure more diffraction-like secondary structures at large \t\ [5, 7, 8]. The chiral-bag 
model with qq condensate cloud, which we studied [9], predicts after the first dip-bump structure 
a smooth approximately exponential fall-off (known as Orear fall-off) and then a slower fall-off due 
to the transition from the nonperturbative regime to the perturbative regime. This change in the 
behavior of dcr/dt was shown only schematically in our previous work [9]. We have now been able 
to quantitatively address this question and study the predicted change of dcr/dt. Results of our 
investigation and the implications for the combined role of perturbative and nonperturbative QCD 
dynamics are briefly reported here. 

We view pp elastic scattering in the perturbative regime as a hard collision of a valence quark 
from one proton with a valence quark from the other proton (Fig.l). The collision carries off 
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Figure 1: Hard collision of valence quarks from two different protons 



the whole momentum transfer. This dynamical picture brings new features in our calculations: 
1 probability amplitude of a quark to have, say, momentum p when the proton has momentum P 
in the cm. frame. 2) Quark-quark elastic amplitude at high energy and large momentum transfer, 
which is in the domain of perturbative QCD. The latter has been the focus of extensive studies 
following the original work of Balitsky, Fadin, Kuraev, and Lipatov (BFKL) [12]. The present 
status is that the qq elastic scattering occurs via Reggeized gluon ladders with rungs of gluons 
which represent gluon emissions in inelastic processes (BFKL Pomeron). It is a crossing-even 
amplitude which is a cut in the angular momentum plane with a fixed branch point at olbfkl = 
1+oj. The value of oj in the next-to-leading order (NLO) lies in the range 0.13-0.18 as argued by 
Brodsky et al.[13]. We refer to the BFKL Pomeron with next to leading order corrections included 
as the QCD "hard Pomeron" . In our investigation, we approximate this hard Pomeron by a fixed 
pole and take the qq scattering in Fig. 1 as 

1 

1*1 t.,, 



f{s,t)=i lm s{se-^Y ^_ 2 , (1) 



where s — (p + fc) 2 , t = —q 2 . The phase in Eq.Q follows from the requirement that T(s, t) is a 
crossing even amplitude. Eq.Q represents the hard Pomeron amplitude in our calculations. If we 
want to describe just asymptotic qq scattering, we have to take into account unitarity corrections 
due to infinite exchanges of this Pomeron. This can be done by taking T(s, t) as the Born amplitude 
in an eikonal formulation [14], which leads to a black-disk description and requires 7 g? > 0. The 
radius of the black disk turns out to be R{s) = r uj Ins. Hence, the parameter in Eq.JQ) has the 
physical significance of a length scale that defines the black-disk radius of asymptotic quark-quark 
scattering. 

We next examine how to obtain the pp elastic scattering amplitude from the process shown in 
Fig. 1. Let s be the square of the cm. energy of the two colliding protons: s = (P + K) 2 . s, of 
course, is the square of the cm. energy of the two colliding quarks. From Fig. 1, we see that 
initially we have a quark of momentum p : \p) with a probability amplitude ipip) in the cm. 
frame in which the proton is moving with momentum P. Similarly, we have a second quark with 
k ) with a probability amplitude (p(k) in the cm. frame in which the other proton 



momentum k 

is moving with momentum K = —P. Thus, the initial state of the two colliding quarks is 

|»)=^|^>V(*)|*) ■ (2) 

After the collision, we have a quark with momentum p '— q ': \p—q) with a probability amplitude 
(pip — q), and a quark with momentum k + q : k + q ) with a probability amplitude <p(k + q). So, 
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the final state is 



\.f) = <P{P- Q) \P~ q)<p{k + q) 



(3) 



The pp elastic scattering amplitude due to quark-quark scattering T qq (s, —q 2 ) from Fig. 1 is then 



T qq (s,-q 2 ) = ^^(p*(p-q)(p*(k + q)/k + q (p-q\f op \p) k)ip(p)ip(k) , 



(4) 



where (k + q 



(p~q\ f op \p) 

depends on the invariants s — (p + q) 2 and t = — <f 2 , we can write 



K) is the qq elastic scattering amplitude. Since this amplitude only 



(5) 



(p-q\ T op \p) kJ=T{s,-q ). 
Eq.(0} then takes the form 

T qq {s,-q 2 ) = ^2^2<p*(p- q)(p(p)f(s,-q 2 ) ip* (k + q)<p(k) 



(G) 



This equation makes it evident that tp*(p — q)ip(p) and ip*(k + q)(p(k) are the nonperturbative "im- 
pact factors" which modify the perturbative qq amplitude T(s,—q 2 ). The right-hand-side (RHS) 
of Eq. © needs to be multiplied by a factor of nine to take into account that there are three quarks 
in each proton[15]. We absorb this factor in the constant 7 gg . 

To see the physical meaning of Eq.©, let us assume that we can approximate qq scattering 
in Fig.l by taking some average value of s : s av - Of course, s av is going to be proportional to s. 
Eq.© then takes the form 



Tqa(s, -q- 



^2<P*{P- q)f(p) T(s av ,-q 2 )^2ip*(k + q)ip(k) , 



(7) 



which shows that the impact factors separate out. Each momentum sum in Eq.(JZJ) can now be 
carried out and yields the form factor associated with the quark probability density in the cm. 
frame. This probability density is Lorentz contracted, which means if pa(r') is the quark probability 
density at r' in the proton rest frame and p(r) is the probability density at r in the cm. frame, 
then 

p(b + e 3 z) = 7p (6 + e 3 7z), (8) 

where 7 is the Lorentz contraction factor: 7 = E/M = y/s/(2M), f=b + e 3 z, and £3 is the unit 
vector in the direction of P, i.e., the z-axis. If F(q) is the form factor associated with po(f): 



F(q) 



d 3 r e^p (r), 



(9) 



and po(r) is spherically symmetric, then 
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(10) 
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In deriving Eci. (|10|) . we have used p(r) = ijj*{r)ijj{f), where the quark wave function tp(r)is related 
to its momentum wave function (p(p) via the plane wave expansion: 

^ = £^($ (ii) 

p 

Eq.Q now takes the form 

T qq ( Sl ~q 2 ) ~ F(q ± ) f(s av , -q 2 ) F(q x ), (® = ^ - 0). (12) 

The structure of Eq. (|T2")) is easy to understand. It is the usual quantum-mechanical scattering 
amplitude of two composite objects described by the form factors and interacting via a basic process 
whose amplitude is T(s av , —q 2 ). We take the form factor F(q) describing the quark probability 
density or number density in the nucleon rest frame to be a dipole: 

m = (i+£*y 2 , as) 

so that it satisfies the dimensional counting behavior t~ 2 for the form factor of a proton made up 
of three quarks [16,17,18]. 

Next we go back to Eq. © and no longer make the approximation of replacing s by an average 
value. Inserting Eq.Q in Eq.JfJJl, we obtain 

T qq (s, -q 2 ) = J2J2 V*(P~ <1)V(P) V (s e-^r * _ 2 f*(k + (f)<p(k) . (14) 

Introducing light-cone variables P+ = Po + P3 , P_ = Pq — P 3 ,p+ = po + P3 ,P- = Po — P3 ,etc. 
and writing p+ = x P+ , fe_ = x 1 K- ,wc find s ~m's, when P+ , K- — ► 00 . Eq. (|14|l then takes 
the separable form 



T qq (s,-q 2 )= |^/(P-?>(P> 1+U | i lm s(se-^r g 2 + r -a ( £ ^ + * + " 



(15) 



In a frame where P+ — > 00 , 

( "l 2, M 2 x 2 1 
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where 



00 

J(g±,o(x)) = /" 6*J (6g±){6iTi[H} 2 - (17) 



Here M is the nucleon mass, mo is the mass parameter that occurs in the form factor Eg. l|13|) . 
1 

a = (^j 1 + M 2 x 2 ^j 2 , and q ~ <fj_ . In deriving Eq. (|16[) . we use momentum wave function (p(p) 
obtained from the Lorentz contracted probability density. It can be related to the rest frame wave 
function <po(j?) in the following way: 

P3 

f(p± + e 3 p 3 ) = <f (p±+e 3 — ), (18) 
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and yields the result 



ip{p ± + e 3 p 3 ) 



2ir ; 



Pi 



Pi 
I 2 



(19) 



(Vo is the quantization volume in the rest frame.) The integral I(q±,a(x)) can be evaluated 
analytically, and we obtain 



I(q±,a 



< a » = 8^ ln( °' + a) + i ln(a ' + a) - ^ 



3a' 3 
ln(a' + a) + — ■ 



(20) 



where a' 2 = 
expressed in the form 



, a 2 = a' 2 + 1 . Let us denote by F(q±) the RHS of Ea.jlfy Eq.||TSJ| can then be 



T qq (s, -q 2 ) = F(q±) i lqq s (s e-*i )" 



1*1+^ 



For a' 2 = >> 1 and a 2 ~ a' 2 , Eq.JU yields 



J(gx,a(i)) 
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(|t| = g 2 ^«i)- 



Substituting this on the RHS of Eq.lJTBJl, we find 



(21) 



(22) 



(23) 



Ea. (|2*T|) then leads to an amplitude 



T qq {s,-q 2 ) 



i "tqqsjs e l ?) u 



This results in differential cross section behavior for fixed s and large \t\: 

1 



da 
~df 



(S » \t\ » 771q + AM 



(24) 



(25) 



Eq. (|25|) shows that we obtain the behavior predicted by the perturbative QCD dimensional count- 
ing rules[16,17,18] for large \t\. 

In our pp elastic scattering model, we now have two hard-collision amplitudes: one due to u> 
exchange, the other due to the hard Pomeron exchange. Both collisions are accompanied by cloud- 
cloud diffraction scattering that reduces these amplitudes by an absorption factor exp(ix(s, 0)) 
[19]. So the sum of the two hard amplitudes becomes 



Tx(s,t) 



± 



F 2 {t) 



^ 2 (gx) 



7S— — -+ij qq s(se _ 
m z - t \t\+ r 



, (+for pp, -forpp). 



Using the same parameterization as before [9] , 



7 e l * {s ' 0) = 70 



Ti 



we find 



Ti(s,t) 



7o + 



7i 
(s e~ % ~ 



± s 



F 2 {t) 
m 2 — t 



+ i % q s (s e 2 



•^ 2 (gx) " 

\t\+To\ 



(26) 

(27) 
(28) 
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where 7gg = Iqq/l- The qq hard scattering term brings four new parameters: i) 7 gg which measures 
the relative strength of this term compared to the w exchange term; ii) oibfkl = 1+w which 
controls the high energy behavior; iii) r which provides the length scale for the black-disk radius 
of qq asymptotic scattering; iv) m which determines the quark wave function ^o{r) = \J po(r) and 
the size of the quark bag. Because of the different physical aspects associated with them, these 
four parameters form a minimal set. 

We determine the parameters of the model by requiring that the model should describe satis- 
factorily the asymptotic behavior of <r to t(s) and p(s) as well as the measured pp elastic dcr/dt at 
Vs=546 GeV [20], 630 GeV [21], and 1.8 TeV [22, 23]. The results of this investigation arc shown 
in Figs. 2-4 together with the experimental data. We obtain quite satisfactory descriptions. The 
dotted curves in Figs. 2 and 3 represent the error bands given by Cudell et. al. (COMPETE 
Collaboration) to their best fit [24]. Wc notice that our a to t(s) curve lies within their error band 
closer to the lower curve, but our p pp (s) curve (dashed curve in Fig. 3) deviates from the band. As 
noted by Cudell et. al., such a deviation is not surprising - since a hard Pomeron occurs in our 
calculations and not in theirs. In fact, this hard Pomeron in conjunction with a crossing-odd ab- 
sorptive correction [19] in our model leads to a crossing-odd amplitude (an odderon) and produces 
a visible difference between pp p (s) and p pp (s) at large y/s. The parameters describing the soft 
(small |i|) diffraction amplitude and the hard (large |i|) w-exchange amplitude have been discussed 
before [9]. Their values are: R = 2.77, Rx = 0.0491, a a = 0.245, ai = 0.126, 770 = 0.0844, c = 
0.00, cr = 2.70, A = 0.727, d Q = 13.0, a = 0.246, % = 1.53, 71 = 0.00, a = 1.46 (the unit of energy 
is 1 GeV). The parameters (3 and m are kept fixed as previously: /3=3.075, to=0.801. There are 
now seventeen adjustable parameters. The four new parameters describing the hard (large \t\) qq 
amplitude have the values ^ qq — 0.03, u — 0.15, r = 2.00, to 2 , = 12.0. (This value of to 2 , leads 
to a valence quark-bag of r.m.s. radius 0.2 F, while that of the baryonic charge core is 0.44 F.) 
These four parameters, however, cannot be determined reliably, because no large |i| elastic data 
are available in the TeV energy region. 

Our prediction for pp elastic differential cross section at LHC at ^fs= 14 TeV for the whole 
momentum transfer range \t\ = - 10 GeV 2 is now given in Fig. 5 (solid curve). We obtain 
for G to t and p pp the values 110 mb and 0.120 respectively. Also given in Fig. 5 are separate 
dcr/dt due to diffraction (dotted curve), due to hard w-exchange (dot-dashed curve), and due to 
hard qq scattering (dashed curve). As expected in our model, we find that in the small \t\ region 
(\t\ ~0 - 0.5 GeV 2 ) diffraction dominates, in the intermediate |i| region (\t\ ~ 1.0 - 4.0 GeV 2 ) u- 
exchange dominates, and in the large \t\ region {\t\ > 6.0 GeV 2 ) qq scattering dominates. The 
three \t\ regions correspond to cloud-cloud interaction, core-core scattering due to w-exchange, and 
valence qq scattering via QCD hard Pomeron. Therefore, they reflect the composite structure of 
the nucleon with an outer cloud, an inner core of topological baryonic charge, and a still smaller 
quark-bag of valence quarks. 

We note that pp elastic differential cross section in the energy range y/s= 27 - 62 GeV and 
\t\ > 3.5 Gev 2 was observed to be approximately energy independent and falling off as i~ 8 . This 
was interpreted as due to the independent exchanges of three perturbative gluons [25, 26]. Later 
it was pointed out that the three gluons would Reggeize, so that color-octet exchanges would be 
suppressed. Instead, three color-singlet exchanges would take their place [27]. Eventually, as \t\ 
increases, a single color-singlet exchange would dominate and lead to a t~ 10 fall-off as predicted 
by the perturbative QCD dimensional counting rules [16,17,18]. In our model, the dimensional 
counting behavior t~ 10 of dcr/dt originates from the hard qq amplitude in Eq. (28). This amplitude 
leads to a distinct change in the slope of the differential cross section from the intermediate \t\ 
region to the large \t\ region as seen in Fig. 5. For example, for 1.0 < \t\ < 3.0 GeV 2 , dcr/dt 
drops by more than two orders of magnitude, while for 7.0 < \t\ < 9.0 GeV 2 , dcr/dt drops by a 
factor of 4.2, i.e. less than an order of magnitude. Similar decrease in dcr/dt slope was observed 
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at ISR by De Kerret et.al. for \t\ > 6.5 GeV 2 at a much lower energy: ^/s= 53 GeV [28]. Lepage 
and Brodsky[18], however, pointed out that at such low energies it would be hard to distinguish 
between amplitudes that lead to i -8 and t~ w asymptotic behavior. 

We conclude that, if precise measurement by the TOTEM group corroborates our predicted 
slow fall-off of pp elastic dcr/dt in the large \t\ region, then that will provide evidence for the hard qq 
amplitude occurring in Eq. (28). This, in turn, will imply: i) presence of the QCD hard Pomeron, 
ii) perturbative QCD dimensional counting behavior at asymptotic \t\ (>> 10 GeV 2 ), and iii) the 
confinement of valence quarks in a small region within the proton. 

Much of this work was done when one of us (MMI) was at the Yang Institute for Theoretical 
Physics at SUNY Stony Brook on sabbatical leave. He wishes to thank George Sterman, Director 
of the Institute, and other colleagues there for their hospitality. He also wishes to thank Michael 
Rijssenbeek and George Sterman for stimulating discussions. 

References 

1. TOTEM: Technical Design Report, Jan. 2004 (CERN-LHCC-2004-002). 

2. A. Donnachie, H.G. Dosch, P.V. Landshoff, and O. Nachtmann, Pomeron Physics and QCD, 
Cambridge University Press (2002). 

3. V. A. Khoze, A. D. Martin, and M.G. Ryskin, Eur. Phys.J C18 (2000) 167. 

4. J. Bartels, E. Gotsman, E. Levin, M. Lublinsky, and U. Maor, Phys Lett.B 556 (2003) 114. 

5. C. Bourrely, J. Soffer, and T.T. Wu, Eur. Phys.J C28 (2003) 97. 

6. H.Cheng and T.T.Wu, Expanding Protons: Scattering at High Energies, MIT Press, Cam- 
bridge, MA (1987). 

7. P. Desgrolard, M. Giffon, E. Matynov, and E. Prcdazzi, Eur. Phys.J C16 (2000) 499. 

8. V. A. Petrov and A.V. Prokudin, Eur. Phys.J C23 (2002) 135. 

9. M. M. Islam, R. J. Luddy, and A.V. Prokudin, Mod. Phys. Lett.A 18 (2003) 743. 

10. M. M. Islam, R. J. Luddy, and A.V. Prokudin, in AIP Conference Proceedings Vol.698, 
2003 (Intersections of Particle and Nuclear Physics, edited by Z. Parsa) p. 142. hep-ph/0307355 

11. M.M. Block, E.M. Gregores, F. Halzen and G. Pancheri, Phys.Rev. D60, (1999) 054024. 

12. J. R. Forshaw and D. A. Ross, Quantum Chromodynamics and the Pomeron, Cambridge 
University Press (1997). 

13. S. J. Brodsky, V. S. Fadin, V.T. Kim, L. N. Lipatov, and G. B. Pivovarov, JETP Lett. 70 
(1999) 155. 

14. A. B. Kaidalov, Regge Poles in QCD in Handbook of QCD, Vol.1, edited by M. Shifman 
and B. Ioffc (World Scientific, 2001) p. 603. 

15. The quarks in our field theory model (Ref. 10) are massless effective color-singlet quarks 
and not the QCD current quarks. 

16. V. A. Matveev, R. M. Muradian, and A. N. Tavkhelidze, Lett.Nuovo Cimento 7 (1973) 719. 

17. S. J. Brodsky and G. R. Farrar, Phys.Rev.Lett.31 (1973) 1153; Phys.Rev.D 11 (1975) 1309. 

18. G.P. Lepage and S.J. Brodsky, Phys.Rev.D 22 (1980) 2157. 

19. M.M. Islam, B. Innocente, T. Fearnley, and G. Sanguinetti, Europhys. Lett. 4 (1987) 189. 

20. M. Bozzo et.al., UA4 Collaboration, Phys.Lett. B147 (1984) 385; B155, (1985) 197. 

21. D. Bernard et.al. Phys.Lett. B171 (1986) 142. 

22. N. Amos et.al., Phys.Lett. B247 (1990) 127. 

23. F. Abe et.al., Phys.Rev. D50 (1994) 5518. 

24. J.R. Cudell et.al., Phys.Rev.Lett. 89 (2002) 201801. 

25. A. Donnachie and P.V. Landshoff, Z. Phys. C 2(1979)55. 

26. A. Donnachie and P.V. Landshoff, Phys. Lett. B 387(1996) 637. 

27. M. G. Sotiropoulos and G. Sterman, Nucl.Phys. B 425 (1994) 489. 

28. H. De Kerret et al., Phys. Lett. B 68(1977) 374. 



7 




- (GeV) 



<i (GeV) 



Figure 2: Solid curve represents our calculated Figure 3: Solid and dashed curves represent our 
total cross section as a function of yfs. Dotted calculated pp p and p pp as functions of y/s. Dot- 
curves represent the error band given by Cudell ted curves represent the error band given by 
et al. [24]. Cudell et al. [24]. 
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Figure 4: Solid curves show our calculated da / dt 
at y/s = 546, 630 and 1800 GeV. Experimental 
data are from references [20], [21] and [22, 23]. 



Figure 5: Solid curve shows our predicted dcr/dt 
for pp elastic scattering at y/s =14 TeV at LHC. 
Dotted curve represents dcr/dt due to diffraction 
only. Similarly, dot-dashed curve and dashed 
curve represent dcr/dt due to hard w-exchange 
and hard qq scattering only. 
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